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Exciton-polariton condensates display a variety of intriguing pattern- forming behaviors, partic- 
ularly when confined in potential traps. It has previously been predicted that triangular lattices 
of vortices of the same sign will form spontaneously as the result of surface instabilities in a har- 
monic trap. However, natural disorder, deviation of the external potential from circular symmetry, 
or higher-order terms modifying the dynamical equations may all have detrimental effects and 
destabilize the circular trajectories of vortices. Here we address these issues, by characterizing the 
robustness of the vortex lattice against disorder and deformations of the trapping potential. Since 
most experiments use time integrated measurements it would be hard to observe directly the rotat- 
ing vortex lattices or distinguish them from vortex-free states. We suggest how these difficulties can 
be overcome and present an experimentally viable interference-imaging scheme that would allow the 
detection of rotating vortex lattices. 
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I. INTRODUCTION 

The experimental realization of Bose-Einstein conden- 
sates of quasiparticles in solid-state systems has pro- 
vided a new example of macroscopically coherent quan- 
tum states of matter. Solid-state condensates exhibit 
striking fundamental differences from traditional quan- 
tum fluids such as atomic Bose gases and superfluid liq- 
uid helium. Much recent interest has been devoted 
to exciton-polaritons, which exist as normal modes of 
strongly coupled excitons and photons in semiconductor 
microcavities. Due to the short lifetime of the quasi- 
particles, these condensates exist in a dynamic balance 
between pumping and decay rather than in true thermal 
equilibrium. As such, these systems are capable of spatial 
and temporal pattern formation similar to that exhibited 
by lasers. -^^'-^^ 

Hallmark features of superfluidity such as quantized 
vortices^'^'^ "'^^'"'^^'^^ and persistent currents^ have been 
observed in polariton condensates. Spatial patterns and 
soliton dynamics have been seen to arise spontaneously in 
both one^ and two^^'^^ dimensions. Intrinsic disorder in 
the material^^ or deliberately designed potentials^ induce 
currents that may further facilitate symmetry breaking 
in the polariton system leading to novel patterns. In 
particular, vortices have been seen to nucleate and get 
trapped because of the nonuniformities of the potential 
landscape. ^^'^^ 

Most of the models that have been successfully and 
extensively used to describe pattern formation in po- 
lariton condensates are based on the complex Ginzburg- 
Landau equation (cGLE), which is also used in the the- 
ory of lasers and thus illustrates the similarities be- 
tween the two systems. Such models have been applied 
to vortex nucleation and motion, ^^'^^^^ spatial density 
modulations, ^^'^^ solitons"*^^ and vortex lattices. ^^^^ 



Various different modifications of this model have been 
considered, such as including separate dynamics of the 
reservoir, or modeling processes of thermal relaxation. 
Under the assumption of fast relaxation of the reser- 
voir and small condensate densities these modifications 
all simplify and the model can be written as: 



(1) 



where V{x,y) describes a trapping potential, a{x,y) is 
a (spatially varying) pump rate, and a is a nonlinearity 
which causes pumping to reduce as density increases — 
some such nonlinear process is clearly necessary for the 
system to be stable, and the form included is the simplest 
such nonlinearity. Linear losses are included in a. An- 
ticipating that we will consider a harmonically confined 
system, Eq. (1) is stated in oscillator units, measuring 
energy in un its of ficj, length in units of the oscillator 
length / = yjhlmuj^ and time in units of cj"^, where 
is the oscillator frequency of the trapping potential, and 
m is the effective mass of the polariton. The values as- 
sumed in the numerical simulations are discussed at the 
start of Sec. III. 

The polarization degree of freedom of the polaritons 
provides additional possibilities in pattern and defect for- 
mation, such as allowing for spatially separate vortices of 
left and right polarization.^ Further possibilities for pat- 
tern formation arise in the presence of a magnetic field 
that favors circular polarization and therefore competes 
with the interactions, which favor equal densities of the 
polarization components, and also with any anisotropy 
resulting from strain fields induced by mechanical stress. 
The cGLE can be modified to include these effects, but 
in this paper we will consider linear polarization only. 
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In the absence of pumping and decay terms, the steady 
states of the cGLE are often weh approximated by the 
stationary Thomas-Fermi solution. One rather surpris- 
ing feature of the cGLE is that there exist many cases 
where the presence of pumping and decay significantly 
modifies the stable density profile. A particularly sur- 
prising case is that in a harmonic trap, the stationary 
state is unstable, and a rotating lattice is formed as the 
result of this instability.^^ Spontaneously formed vortices 
have been observed in polariton condensates, ^'^'^^"^^ but 
these have been isolated, stationary vortices, pinned by 
the disorder potential in the sample. The spontaneous 
formation of a regular, rotating vortex lattice^^'^^ has 
not been seen experimentally. The aim of this paper is 
therefore to study the robustness of such lattices, and to 
discuss how one may overcome the difficulties that might 
arise in detecting a rotating vortex lattice. 

We consider the robustness of vortex lattices with 
respect both to deformations of the circular harmonic 
trapping potential considered previously,^^'^^ and also to 
modifications of the cGLE, such as energy relaxation. 
We find that there is definite sensitivity to the trap- 
ping potential — both elliptical distortion as well as back- 
ground disorder can destroy the regular vortex lattice and 
replace it with a chaotic regime, and we present numer- 
ical results for the critical distortion and disorder that 
destroy the lattice. This fragility with respect to the ge- 
ometry of the trapping potential contrasts with the rela- 
tive robustness of the vortex lattice to other terms that 
might be present in the cGLE, such as phenomenological 
approaches to the inclusion of relaxation processes. 

Detection of moving vortices in a polariton conden- 
sate is non- trivial. Polariton condensates can be imaged 
directly using the photons escaping from the semicon- 
ductor microcavity. Each image requires that light be 
gathered over a non- negligible period of time, and there- 
fore any moving feature will tend to smear out. Hence 
rotating configurations would not be directly visible in 
time integrated images. We show that these problems 
can be overcome by using a defocused homo dyne imag- 
ing scheme, allowing stationary repeat able images to be 
observed, which fully characterize the vortex lattice. 

The remainder of this paper is organized as follows. 
Section II discusses the instability of the stationary state, 
showing how this instability can be found analytically 
for weak pumping and decay. Section III then considers 
the vortex lattices that result from this linear instabil- 
ity, presenting results for elliptical trapping potentials, 
and traps with additional disorder present. Section IV 
presents schemes for detecting rotating vortex lattices, 
and section V summarizes these results. 



II. INSTABILITY OF THE NON-ROTATING 
STATE 

As commented above, the presence of pumping and de- 
cay destabilizes the stationary Thomas-Fermi potential 



in a circular trap, which in turn leads to the appearance 
of the vortex lattice. In this section we discuss how, for 
weak pumping and decay, this instability can be found 
analytically, providing some insight into how and why 
the stationary state becomes unstable. This section pro- 
vides an outline and the main results of the analytical 
stability calculation; further technical details are given 
in Appendix A. 

A. Linear stability analysis 

In order to find the leading order effects of pumping 
and decay, one may linearize both in fluctuations and 
in the effects of pumping and decay, i.e., one considers 
first the density profile and fluctuations about this in the 
absence of pumping and decay, and then determines how 
these are modified by non-zero a and a. 

Using %l) = y/pe^^ and V(x, y) = x'^ -\- y'^ = r'^ , one may 
rewrite Eq. (1) as a pair of real equations: 

9tP + V • (pv) = {a- ap)p 
2at^+V(p + r2 + |^P) =0, 

where v = Vcj). In these equations, an approximation has 
already been made in neglecting the quantum pressure 
terms; this approximation is appropriate for sufficiently 
smooth density profiles (i.e., for clouds much larger than 
the healing length = l/^/p). By substituting p 
p -h h{t),v ^ V -\- u(t) into Eqs. (2), and linearizing in 
1?, /i, the fluctuations obey 

dth = —V • {pu + hv) -\- {a — 2a p) h 
dtu = -V{h/2 + V'u). 

In the absence of pumping and decay, the steady state 
is given by p = ptf = (/i — r^)6(/i — r^), where 6 is 
the unit step function and = 0, so the terms in Eq. (3) 
involving v vanish. In this case Eq. (3) can be reduced 
to the single equation 

d?h _ I d ( dh\ P . 

^ " 2^^ V ; ^ 2^d0^ ' ^ ^ 

which has eigenf unctions idthn,s = ^n,shn,s labeled 
by angular momentum hn,s{r,0) = e*^*^/in,s(r), where 
hn,s{'^) are hypergeometric functions of the radial vari- 
able r (see Appendix A 2). These eigenfunctions have 
frequencies 

C0n,s = \/5(l + 2n) + 2n(n-hl). (5) 

B. Perturbative shift with pumping and decay 

We may now use the solutions (4) and (5) as a basis 
on which to consider perturbatively the effects of pump- 
ing and decay. We will work here to first order in a and 
(7, in order to see the qualitative change introduced by 
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pumping. This first-order treatment will predict an in- 
stability that already explains the behavior seen when all 
orders are considered in the numerical results presented 
in Ref. 26 and in Sec. III. 

At first order, the steady state of Eq. (2) is given 
hy p = ptf and v = —apTFr/6. Integrating the first 
of Eq. (2) over the entire space gives p = 3a/2a. To 
find the change to the frequencies of the fluctuations, it 
is helpful to consider Eq. (3) in the form 

{Looi + ScoiXi^^i + 5rpf) = i{C + a<5/:)(V'o1 + Hf\ (6) 

where ijjf' represent the right eigenfunctions of the matrix 
of differential operators iC. Note that in this identifica- 
tion, the terms involving v in Eq. (3) should be considered 
part of The change to the eigenfrequency induced by 
pumping and decay is thus 



Suji = ia 



I 



J 



(7) 



where ipf indicate left eigenfunctions. Finding left and 
right eigenfunctions separately is necessary as the oper- 
ator iC is not self-adjoint; further details are given in 
Appendix A 1 . 



1. Instability of the stationary solution 

Using the hypergeometric solutions described above, 
and following the outlined procedure (see Appendix A 3 
for details), one finds the shift 



s2-s{l + 2n)-2n(n+l) 



+ 5(1 + 2n) + 2n(n + 1) 



(8) 



The crucial feature of this result is that for any n there 
is always s (for instance s ^ n), such that this expres- 
sion is positive. This implies that the eigenfrequencies 
acquire a positive imaginary part and the corresponding 
eigenmodes, which are those with large angular momen- 
tum, will therefore grow, leading to vortex nucleation 
(cf. atomic Bose-Einstein condensates where vortices are 
nucleated when external stirring causes unstable eigen- 
modes to arise^^). 



states. At lowest order of approximation this leads to 
the same ansatz.^^'^^ Such terms have also been used 
in modeling polariton condensates. ^^'^^'^^ As described 
in Appendix A 4, including such an effect at linear or- 
der changes the steady state chemical potential to p = 
3ce/(2cr -h 377), and the frequency shift becomes 



5ui 



2cr 

2cr+3?7 



-5(l + 2n) -2n(n + l) 



-h s(l -h 2n) -h 2n(n -h 1) 



(9) 



Thus such relaxation, although driving the cGLE to- 
ward its low-energy states, is not sufficient to kill the 
instability — for any n and any finite r^/a, there always 
exist modes (for instance, with s ^ ncj/r]) which will 
become unstable. 



III. PROPERTIES OF THE ROTATING STATE 

The above stability analysis assumed an infinite pump 
spot, however the same instability still exists as long as 
the pumping spot is large r than the condensate cloud size 
(given by ttf = ^/p = \/3ce/2cr) . As found previously in 
Ref. 26, the instability leads to a rotating vortex lattice 
in a circular pump spot. In this section we discuss the 
appearance of rotating states in cases where the trapping 
and pumping deviate from the ideal circular harmonic 
case considered in previous works. ^^'^^ To discuss such 
cases, we present results of numerical simulations. 

In choosing values of the parameters a and a we follow 
the numbers used in Refs. 26 and 28. The parameter a 
may be directly related to the blueshift of the condensate 
when pumping at a given multiple of the threshold power. 
In this way a is estimated in Ref. 26 as cr = 0.3. The di- 
mensionless gain a is given by the background linewidth 

and the pumping power 7 as a = — K)/huj. Esti- 
mating the polariton lifetime at r 5 ps gives a linear 
decay rate = h/r = 0.13 meV. To facilitate comparison 
with Refs. 26 and 28, we choose a = 4.4. In the stress 
trap constructed by Balili et al,^''^^ huj = 0.066 meV for 
the harmonic confinement under the stressor, and the ef- 
fective mass of the polariton is m = 7 x lO^^mg. For this 
value of huj our chosen a roughly corresponds to pumping 
at twice the threshold. 



2. Including phenomenological relaxation 

The same procedure may be repeated includinging var- 
ious modifications to the cGLE. One such modification 
that may be easily incorporated is an energy-relaxation 
term, causing time evolution toward lower energy states, 
i.e., dt^p ^ (1 + ir])dtip. The introduction of this term 
is common in literature on atomic Bose-Einstein conden- 
sates, where it models the atomic transfer between the 
condensate and the thermal cloud. These effects can 
be accurately modeled by a quantum-Boltzmann equa- 
tion describing the population dynamics of the quantum 



A. Non-circular geometry 

Earlier considerations of spontaneously formed rotat- 
ing vortex lattices in harmonic potentials'^ '^^ have as- 
sumed circular symmetry. Here we lift this assumption 
by considering an elliptical geometry. For simplicity in 
illustrating the effects of elliptical geometry, we take both 
the trapping potential and the pumping spot in Eq. (1) 
to be described by a parameter ^, such that^^ 

V{x,y) = ^ + SV- (10) 



FIG. 1. (color online) Density (grayscale) and superfluid 
streamlines (arrows) for different sizes of the pumping spot 
(red ellipse). Left: {S = 0.9, R = 7) The surface instability 
leads to a spontaneously formed lattice with vortices mov- 
ing on elliptical trajectories. Right: {6 = 0.9, = 4) The 
pumping spot is smaller than rxF and the instability does 
not appear. 



Using parameter values as detailed above, the pumping is 
taken to be constant inside the ellipse /S'^ -\-S'^y'^ = R^. 
The finite lifetime of the polaritons leads to a linear loss 
outside the pumping region. Since we are considering 
pumping at twice threshold, 7 = 2/^, the value of a out- 
side the pump spot is the same as inside but with oppo- 
site sign. As in the circular case,^^ the surface instability 
will appear only if the pumping spot is sufficiently large, 
as illustrated in Fig. 1 for ^ = 0.9. For our subsequent 
analysis, we keep a fixed R = 7. 

The movement of individual vortices in the lattice can 
be traced using a minimum-finding algorithm. The orbit 
is obtained by averaging the trajectory over a number of 
revolutions. In the regular vortex lattice, the individual 
vortices move on elliptic orbits as shown in Fig. 2. Be- 
cause the vortices have a definite sense of rotation, the 
motion is not symmetric around the major axis, and the 
resulting elliptic vortex lattice is tilted with respect to 
the trap. 

In a sufficiently eccentric geometry, a steady array of 
vortices fails to form. We find a critical 6 ^ 0.85 for which 
some vortices move on elliptical orbits [Fig. 3, left]. At 
S = 0.8, the behavior is entirely chaotic, with vortices 
continually entering and leaving the cloud [Fig. 3, right], 
causing fragmentation of the condensate. 

In the limit of extreme eccentricity, the condensate 
tries to establish a density profile of approximately uni- 
form width in the transverse direction, leading to a less 
fragmented condensate at (5 < 0.2 [Fig. 4]. The instability 
towards entrance of vortices is still present where the spot 
is widest: vortices entering and leaving the cloud cause 
transverse oscillations. As illustrated in Fig. 4, vortices 
of opposite signs may be present. The distribution of 
vortices of different sign may lead to a net effective pres- 
sure in the condensate pushing the perturbation along 
the length of the condensate until the pressure gradient 



FIG. 2. (color online) Left: the averaged orbits of two vortices 
(green, dashed and magenta, dash-dotted lines, respectively). 
The red, solid ellipse marks the pumping spot. Right: snap- 
shot of the vortex lattice with the tracked vortices marked by 
green and magenta crosses. 




FIG. 3. (color online) Left: S = 0.85; the regular vortex lattice 
is broken up, with some vortices still on regular elliptic orbits. 
Right: 6 = 0.80; as 6 decreases, the rotating lattice gives way 
to a fragmented state with vortices continually entering and 
leaving the cloud. 



disappears. This effect is illustrated in Fig. 4 (bottom), 
showing the final position of the region of perturbations. 
Note that a true stationary state is not obtained: oscil- 
lations are still present. 

We see that the formation of a regular vortex lattice 
is sensitive to the geometry of the system, in contrast to 
the robustness against refinements of the pumping model 
discussed in Sec. II B 2 and Sec. IIIC. While a rotating 
vortex lattice will still form under small deviations from 
circular geometry, moderate eccentricity will prevent the 
lattice from forming. 



B. Effects of disorder on the rotating state 



In the samples used in experiments disorder is always 
present. Here we take this into account by including a 
disorder term in V{x,y). This is modeled as a Gaussian 




FIG. 4. (color online) Top left: 6 = 0.25. In a highly ec- 
centric trap, the density tries to establish a profile of approx- 
imately uniform width along the length of the trap. Insta- 
bility towards entrance of vortices leads to transverse oscilla- 
tions. Note that the vertical axis has been rescaled so that 
the pumping spot appears circular in order to accommodate 
the extreme eccentricity. Top right: 6 = 0.1. Extreme ec- 
centricity suppresses fragmentation of the cloud. Transverse 
oscillations develop at the midpoint. Bottom left: Sponta- 
neously entering vortices cause an effective pressure gradient, 
pushing the region of perturbations towards one end of the 
trap. The final position is shown. Bottom right: Vortices 
of both signs are present in the perturbed region. Arrows 
indicate circulation. 



random field with real-space correlation function 



(^dis(r)Fdis(rO)=^o'exp - 



2^ 



(11) 



where the amplitude is a measure of the strength of 
the disorder, and ^ is the correlation length. (Such a po- 
tential can be constructed by sampling the Fourier com- 
ponents of the disorder potential from a Gaussian distri- 
bution with a /c-dependent variance; see Appendix B.) 

Varying Vq and ^ we find that the formation of a vor- 
tex lattice is sensitive to the amount of disorder, with the 
lattice being destroyed for Vq > Krit — 6.5. When the 
disorder strength exceeds 14rit the behavior is chaotic and 
no stationary state forms. However, the critical amount 
of disorder is only weakly dependent on the disorder cor- 
relation length ^, as illustrated by the nearly horizontal 
border in the (l^, ^)-diagram of Fig. 5. For the stress 
trap of Ref. 2 our estimate translates to a critical disor- 
der Vcrit — 0.4 meV, subject to uncertainties in estimates 
of experimental quantities as well as effects not accounted 
for in our model, such as relaxation mechanisms. 

It is worth noting that V^rit > ^ give a false 
impression of a very large critical disorder strength. The 
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FIG. 5. (color online) The rotating vortex lattice forms also in 
the presence of weak disorder (a), but is destroyed by a disor- 
der stronger than Vcrit, leading to a chaotic state (b). Bottom: 
the boundary between vortex-lattice and chaotic regimes in 
the (Vb, parameter space. The choices for (Vb, leading to 
states (a) and (b) are indicated in the diagram. The vortex 
lattice is strongly sensitive to the amount of disorder, whereas 
the correlation length is of only marginal importance. 



relevant energy scale is given by effective depth of the 
harmonic potential over the extent of the cloud, and this, 
in turn, is determined by the chemical potential hw. 



C. Effects of other modifications to the cGLE 
model 

There are a number of other possible modifications to 
the simple cGLE model (1) that could potentially desta- 
bilize a vortex lattice or prevent its formation: pres- 
ence of repulsive interactions between condensate parti- 
cles and reservoir particles, presence of higher-order non- 
linear terms, or a finite gain linewidth leading to the ap- 
pearance of superdiffusion in the governing equations. In 
this section we comment briefly on possible consequences 
of considering these effects in the cGLE model. 

Repulsive interactions with an excitonic reservoir. Re- 
pulsive interactions between reservoir polaritons and the 
condensate have been predicted theoretically^^ and de- 
tected experimentally.^^ For our geometry, such interac- 
tions may be modeled by adding a Gaussian potential 
to the center of the trap potential V{x^y). The pres- 
ence of this additional potential term does not prevent 
lattice formation, but increases the probability of having 
a central, stationary vortex, as the reservoir provides a 
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pinning site. The presence of Gaussian reservoir repul- 
sion may further change the structure of the lattice: for 
example a sharp and narrow Gaussian favors a square, 
rather than triangular, lattice. 

Higher order nonlinearities. One may consider higher- 
order corrections to polariton interactions in the form 
of a quintic nonlinearity in the cGLE. In some cases, 
inclusion of quintic terms in cGLE models is known to 
change the stability of the solutions. We have checked 
that the introduction of such terms in our model does not 
have any pronounced effect on vortex lattice formation. 

Superdiffusion. Lasers emit particular transverse 
modes that depend on the detuning between the longitu- 
dinal modes of the reservoir (controlled by the resonator 
length) and the frequency at which gain is maximum. -"^^ 
If such gain selection is relevant to the polariton con- 
densates the right-hand side of Eq. (1) acquires a term 
iz/(V^ + A)^?/;, where A is the detuning of the gain above 
the lowest transverse mode. Such mechanisms have not 
been extensively discussed in the context of polariton 
condensation, but the essence appears for example in 
Ref. 40. Our numerical simulations show that although 
the lattice survives for small values of u the rings of vor- 
tices are shifted to the boundary of the condensate. The 
lattice disappears for u > 0.1. 




FIG. 6. Example of a vortex-lattice solution to Eq. (1). Top 
left: density, showing the vortex lattice. Top right: the cor- 
responding spectrum. Note the three bands corresponding to 
the regions of the condensate separated by the rings of vor- 
tices. Bottom: complex phase of the wavefunction (left), and 
the time derivative of the phase (right). Note the separa- 
tion into regions with approximately constant time derivative 
corresponding to the bands in the spectrum. 



IV. DETECTING ROTATING VORTEX 
LATTICES 

We now discuss how time-integrated imaging tech- 
niques can be used to observe a non-stationary vortex 
lattice in experiments. Images are produced by gather- 
ing light escaping from the condensate over a period of 
time long enough that the movement of the vortices is 
non- negligible. We will construct a way of extracting an 
image that remains stationary during the measurement 
despite this movement. In practice, imaging is done by 
combining images from several realizations of the exper- 
iment. It is therefore necessary that the same image is 
produced in each realization. In the following, we discuss 
the relative merits of heterodyne and homodyne imag- 
ing techniques, before combining the desirable features 
of both into an experimentally viable scheme described 
in Sec. IV C. 

In the imaging schemes, we exploit the fact that a 
rotating vortex lattice generates sharp sidebands in the 
spectrum. A representative example is shown in Fig. 6. 
Each sideband corresponds exactly to one of the regions 
of the condensate separated by the rings of vortices. 
By interfering each region with light with matching fre- 
quency, the phase difference between light from the con- 
densate and the interference beam will be constant in 
time, creating a fixed interference pattern, which can be 
integrated over time without smearing. (See Appendix C 
for technical details.) 



A. Heterodyne imaging 

We first consider a heterodyne interference scheme, in 
which the image of the condensate is interfered with a 
monochromatic reference beam without spatial structure, 
i/^ref = ^0 exp[— i(^o + '^0^)]- If '^o is chosen to match one 
of the sidebands in the spectrum, the phase difference 
between reference beam and condensate is independent of 
time in the region corresponding to the chosen sideband. 
Thus regions of constructive and destructive interference 
in this part of the image remain stationary. 

Fig. 7 shows interference images obtained from the ex- 
ample state of Fig. 6. The vortex lattice appears in the 
image as interference fringes, the number of fringes corre- 
sponding to the number of vortices inside the stationary 
image. 

This heterodyne scheme serves as a conceptually 
straightforward illustration of the type of image we would 
like to obtain. From an experimental point of view, 
however, the scheme is impracticable. A fundamental 
problem is that the initial phase is arbitrary. This 
means that while the interference fringes remain station- 
ary within one measurement, there is no mechanism for 
ensuring that they appear in the same place in any other 
realization (i.e., for locking the initial phase difference). 
This is crucial, as it means that while each image, in prin- 
ciple, will be clear, averaging over several realizations will 
cause the interference fringes to disappear. In addition to 
this fundamental problem, there may be practical diffi- 
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FIG. 7. Time- integrated interference image using the het- 
erodyne scheme and the vortex-lattice from Fig. 6. Left: 
uo = 26.5. Right: luq = 14.5. 




FIG. 8. Snapshot of intensity (left) and phase (right) after 
filtering from Fig. 6 to the band 26.25 < Ct; < 26.75 



culties in establishing a sufficiently stable phase reference 
at exactly the required frequency. 



B. Frequency- filtered homodyne imaging 

Since the heterodyne scheme has problems with the 
arbitrary phase reference, we next consider using a ho- 
modyne scheme to avoid this. The simplest homodyne 
scheme that would reveal information about a vortex lat- 
tice is one where an image of the condensate is filtered 
to a narrow frequency band and then interfered with a 
retroreflected and displaced copy of itself. This is a con- 
ceptually straightforward extension of well-established 
experimental techniques for detecting vortices in polari- 
ton condensates.^ 

Frequency filtering [Eq. (C3)] ensures that even after 
displacement, each part of the image is interfered only 
with light of the same frequency, which is necessary to 
create a stationary interference pattern. Density and 
phase of the example lattice from Fig. 6 after filtering 
to 26.25 <uj < 26.75 are shown in Fig. 8. 

The interference pattern created from Fig. 8 depends 
on the amount and direction of the displacement between 
the copies [Fig. 9]. The images can be understood by 
considering how overlapping areas of the image in Fig. 8 
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FIG. 9. Frequency-filtered homodyne scheme applied to the 
solution in Fig. 6 (the result of the filtering shown in Fig 8). 
From top left to bottom right, the displacement [see Eq. (C2)] 
is Xd = 0.2, Xd = 0.5, Xd = 1.0, and Xd = 2.0. {yo = 0.) 



and its retroreflected copy move in and out of phase as 
displacement is varied. In Fig. 9, this leads to destructive 
interference for small displacement, and to a complicated 
interference pattern as displacement is increased. 

Since the image is produced by interfering light from 
the condensate with itself, there is no arbitrary phase 
difference. The resulting image is therefore identical be- 
tween realizations, and images from different realizations 
can be added without smearing. However, interpretation 
of the images is difficult. It is not straightforward to infer 
the vortex lattice from the interference image, let alone 
glean more detailed information about its structure. 



C. Defocused homodyne imaging 

We resolve the problems identified above by consider- 
ing a defocused homodyne scheme constructed to com- 
bine the desirable features of the heterodyne and ho- 
modyne schemes. This will use a spatially uniform and 
nearly monochromatic reference beam, but with a mech- 
anism for locking the arbitrary phase reference between 
realizations. 

From the image of the condensate, light is taken from 
a reference point ro and then defocused and interfered 
with the full image. Here we assume a reference point 
without spatial extent, corresponding to one grid point 
in the numerical simulation. In practice, the reference 
point would be a fixed spot small enough that the spatial 
variation of the phase across it is negligible. Optionally 
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FIG. 10. (color online) Top left: Snapshot of density with 
crosses showing the points from which the reference is taken. 
Top right: Interference with the reference beam taken from 
the red cross, mapping out the interference fringes of all 18 
vortices enclosed by that part of the condensate. Bottom 
left: Interference with light taken from the blue cross, giving 
six stationary interference fringes. Bottom right: Same as 
bottom left, but with the image filtered to 14.15 < cj < 14.75 
before the reference is extracted and allowed to interfere with 
the filtered image. 

a frequency filter [Eq. (C3)] may be applied as the first 
step of the scheme. 

The scheme then amounts to forming the interference 
/(r, t) = |'^(r, + ilj{fo^t)f. Choosing fo to be in a ring 
where the time derivative of the phase is approximately 
constant (see Fig. 6) creates a reference beam with the 
frequency of a side band in the spectrum. Therefore, the 
interference image is qualitatively the same as that pro- 
duced by the heterodyne scheme, but since the phase is 
now provided by the condensate itself, is no longer ar- 
bitrary. By construction, the interference at fo is always 
constructive, thus locking the position of the interference 
fringes. 

Fig. 10 shows the defocused homodyne scheme applied 
to Fig. 6. The reference point fo is indicated on the snap- 
shot of the density. Selecting fo between the vortex rings 
yields an interference image with six fringes correspond- 
ing to the inner ring of six vortices. Selecting fo near the 
edge of the cloud reveals all 18 vortices. 

The bottom panels of Fig. 10 show the effect of fre- 
quency filtering prior to applying the defocused homo- 
dyne scheme. Without filtering, interference fringes are 
mapped out the same way as in the heterodyne scheme. 
However, filtering may yield a clearer picture with better 
contrast. 



FIG. 11. Interference with light taken from the point indi- 
cated by the red cross in Fig. 10 after filtering to a narrow 
frequency band. The panels show the resulting interference 
image for different in-plane momenta in the x-direction. From 
left to right /c^ = 1.0, 2.0 and 5.0. 




FIG. 12. (color online) Top and bottom left: Density and the 
corresponding spectrum. Crosses mark reference points for 
the filtered defocused homodyne scheme. Remaining panels 
show interference images corresponding to light taken from 
the red (top) and blue (bottom) crosses. From left to right 
in both rows, in-plane momentum in the x-direction kx — 0, 
kx = I and kx — 2. 

Slight misalignment (intentional^ or accidental) in the 
interference optics may introduce a phase gradient in one 
of the pathways. In order to investigate whether the 
proposed scheme is robust against such misalignment we 
introduce a finite in-plane momentum k in the reference 
beam, such that '0ref(^, ^) = '^(fo, ^) exp(— • f). This 
results in a distortion of the interference fringes [Fig. 11]. 
The distortion is small as long as k is not too large. 

The distortion is more severe for smaller vortex lat- 
tices, as demonstrated in Fig. 12. Note that this lattice 
rotates in the opposite direction compared with Fig. 6, 
and that this is reflected in the distortion of the fringes 
in the presence of the phase gradient (cf. Fig. 11). 



V. CONCLUSIONS 

We have shown how the instability of high-angular mo- 
mentum modes leading to spontaneous formation of vor- 
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tex lattices can be understood analytically. We have fur- 
ther addressed several difficulties involved in producing 
and detecting vortex lattices. 

Our analytical treatment shows that the surface-mode 
instability is robust against including relaxation effects in 
the simple models of pumping previously considered. Nu- 
merical analysis of the complex Ginsburg-Landau model 
reveals that both non-circular geometry and disorder may 
prevent formation of a vortex lattice, in contrast to its 
robustness against other modifications, such as relax- 
ation, interactions with reservoir polaritons and higher- 
order nonlinear it ies. Sufficient elliptic deformation or 
disorder results in a chaotic, nonstationary state. How- 
ever, small deviations from circular symmetry as well 
as a small amount of disorder do allow a regular vor- 
tex lattice to form. The critical strength of disorder for 
a = 4.4, cr = 0.3 is estimated at ~ G.Sficj, which corre- 
sponds to ~ 0.4 meV at twice- threshold pumping in the 
stress trap of Ref. 2. 

We have further shown that spontaneously formed, ro- 
tating vortex lattices in a polariton condensate would be 
possible to detect using an experimentally feasible defo- 
cused homodyne scheme to generate stable interference 
images. By using a defocused homodyne scheme one can 
ensure that the same image is generated in each real- 
ization of the vortex lattice. Furthermore, the images 
contain direct information relating to the structure of 
the vortex lattice, allowing one in principle to establish 
the number of vortices and their arrangement within the 
lattice. 
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Appendix A: Linear stability analysis 

This appendix provides further details of the analytic 
stability analysis presented in Sec. II. Sec. A 1 discusses 
the perturbation theory for Eq. (3), while Sec. A 2 gives 
the exact zeroth order solutions of Eq. (4). Sec. A3 then 
shows how the resulting integrals can be analytically eval- 
uated, and Sec. A 4 discusses the effects of relaxation 
terms that may be added to the cGLE. 



1. Perturbation theory for Eq. (3) 

As briefly noted in Sec. II, the zeroth order differen- 
tial operator iC is not self adjoint, and so evaluation of 
perturbation theory requires one to separately determine 



the left and right eigenfunctions of this operator. With 
such eigenfunctions identified, one may write the linear 
part of Eq. (6) as 

Following standard perturbation theory, one then writes 



+ j d\^^,{iC-uj^i)5^t (Al) 



The last term will cancel as long as the operator iC has 
left eigenfunctions which satisfy 



1 = j d'^rxP^^iCiP = uJoi j d^ 



for all ijj. 

The right eigenfunctions are given directly by solutions 
of Eq. (3). To find the corresponding left eigenfunc- 
tions, it is convenient to write out the components of 
the eigenfunctions as t/^^ = {hQi,u^Q-,UQ'^.)^ and i/jq- = 
(/iq^, i^^oi' ^OoJ where uq are radial and polar compo- 
nents of the velocity vector u. The zeroth order right 
eigenfunctions obey Eq. (3) with = 0, = 0, cr = 0. 
This can be written out in components as 



R _ i dhoi , , r,,R - ^ 

Writing out the components = (h^Ur^ue)^ , one has 
i d 



(A2) 



r dr 



i dh 
2d^ 



' Oi 



2r 



Integrating by parts (noting that the integral range is < 
r < -y/JI and that rp{r) vanishes at the boundary), one 
may show that if {h^i^u^^^.^u^J = {h^^ ,2pu^;^,,2puf^;) 
then 

I = u;oiJ d^r [h^;h + 2p {u^;^,Ur + ^^o,^^)] 

as required. Hence the above gives the relation between 
the left and right eigenfunctions. 



2. Exact zeroth order solutions 

With the left eigenfunctions as defined above, the shift 
in frequency can be explicitly determined by evaluating 
Eq. (Al). The functions u^q-^u^^. are given in terms of 
Hq- in Eq. (A2), and substituting into the equation for h^- 
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yields Eq. (4). For p = pxF, and hg{r,0)^ = e''^hg{r), 
this equation is hypergeometric. By restricting to solu- 
tions that remain bounded as /i, one finds the solu- 
tions to be the finite order hypergeometric polynomials, 
with frequencies as given in Eq. (5) and wavefunctions 

h^.{r, 0) = e''^r'2Fi (-n, 5 + 1; n + 5 + 1; r^) . (A3) 

This simple form is the two-dimensional analogue of the 
results given in three dimensions in Refs. 42 and 43. Note 
that the two-dimensional harmonic trap is not equiva- 
lent to the highly anisotropic pancake limit of the three- 
dimensional trap — the two dimensional case presented 
here has a hard wall boundary condition in the z di- 
rection. As a result, the excitation spectrum given here 
is not the same as in Ref. 43. 



3. Evaluating overlaps for hypergeometric 
functions 



Using the results above, one finds that 



r ar 



2u;^ dr dr \ dr J 2uj'^r'^ dr 



(A4) 



By noting that v{r) = \v{r)\ vanishes at the limits, one 
may integrate the second and third terms by parts. Using 
Eq. (4), one then finds that Eq. (A4) simplifies to 

J <fr^^,5C^^, = j <fr[a-2a{^Ji-r^)]\h^^\ (A5) 

Hence, making use of /x = 3a/2(7, Eq. (7) becomes 5uJi = 
(ice/2) {S{r'^)i/iJ. - 2) , where 

This integral can be found straightforwardly by using the 
identity 



4. Stability of the generalized order parameter 
equation 

Including the phenomenological relaxation term dis- 
cussed in Sec. JIB 2 the cGLE becomes 



,dij 1 



-v^ + y(x,2/) + |^p 



(A6) 



and the Madelung representation, Eq. (2), becomes 

9tp + V • (pv) = {a-ap^ 2r]dt(t))p 
25t0+(p + r2 + H')=O. ^ ^ 

The second of these equations has now been written ex- 
plicitly for the phase, rather than for v = V0, since 
the phase now appears explicitly in the first equation. 
In keeping with the hydrodynamic approximation which 
neglects V^y^/y^, a similar term rjdtp/p has also been 
neglected. 

Since the second of Eq. (A7) is unchanged, the station- 
ary density profile remains as before {v is of first order in 
pumping and decay terms), with p = —29^0. The con- 
tinuity equation is however modified, so that integrating 
it over the entire space gives / 27rrdr[a — a{p — r'^) — 
r]p]{p — r^) = 0, yielding p = 3a/ {2a + 3?^). 

The linearized fluctuation equations are complicated 
by the explicit appearance of so it is necessary to ex- 
pand (j) ^ (j) -\- ip{t)] one finds 



+ 



h\_(-W- {pu) 
V • {hv) -\- {a — 2a p — r]p) h 



-V • u 



(A8) 



where u = Vif. In the above, the time derivatives have 
been kept together on the left hand side. One may then 
proceed by multiplying both sides by the inverse of the 
matrix appearing on the left, and expanding to first order 
in the pumping and decay terms a, a, and rj. The net 
result of this procedure is that the expression {a — 2a p), 
which appeared in Eq. (A5), is now replaced by {a — 
2a p) ^ {a — 2a p — r]p — rjp). Hence the frequency shift 
is given by SuJi = {i/2) [a - 2{a + r])p + {2a + r]){r'^)i] . 
Combining this with the values of p and {r'^)i/p given 
above gives Eq. (9). 



dh 
pr—- 
dr 



^ P 2 2 

uj pr 



(making use of Eq. (4)). Evaluating the second term, 
and using the eigenvalues equation, one finds that (r^)^ = 
/i(5^+cjf )/(5^ + 2a;|). Using Eq. (5) for the eigenfrequen- 
cies, one recovers the result in Eq. (8). 



Appendix B: Disorder potential 

This appendix gives technical details for the inclu- 
sion of disorder discussed in Sec. IIIB. The disor- 
der is included in the external potential by letting 
V{r) -^V{r)-\- Vdis(^), where Vdis(^) is a Gaussian ran- 
dom field, constructed on the computational grid as a 
Fourier sum Kiis = SzjKj/^) exp (i {kiX -h kjy -h (t)ij)), 
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whose amplitudes a^j = a_i_j are sampled from a Gaus- 
sian distribution with mean zero and standard deviation 



Sij = vo exp [-ll {kj -h . 



(Bl) 



The phase is a random number in the interval [0, 2tt\. 
The quantities V[) and /c, defined for computational con- 
venience, are related to the physical disorder strength Vq 
and correlation length ^ by noting that 



(Kiis(rlKiis{r')) 



and identifying ^ 
where A/c is the 



4(AA:)2 2/2 



exp 



r') 



(B2) 



Eq. (B2) follows from the fact that {a. 



placing /c-space sums by integrals. 



= 2lc and Vo = VV2/(2/cAA:), 
numerical Fourier mode spacing. 

2 and re- 



Appendix C: Constructing interference images 

This appendix gives additional technical detail on the 
calculation of interference images. Vortex-lattice solu- 
tions are generated by propagating Eq. (1). The gen- 
eral principle for interference imaging of rotating vortex 
lattices discussed in Sec. IV does not, however, depend 
on the specific Ginzburg-Landau model. The disper- 
sion spectrum of the condensate is found as the modulus 



squared of the 3D Fourier transform of the condensate 
wavefunction sampled over a period of time: 



(CI) 



The time-integrated interference images are found as 
the integral 



/(r) = J dt |V^(f-hfd,t) 



■'0ref(^-^d,^)r 



(C2) 



where the time T is the time over which light is gathered 
in the measurement. Here we choose T to correspond to 
a few revolutions of the vortex lattice. Since the interfer- 
ence image is stable, it is insensitive to the choice of T. 
The form of V^ref is given by the choice of imaging method 
[Sees. IVA-IVC]. The displacement fd is non-zero only 
for the homodyne scheme in Sec. IV B. 

For the frequency-filtered homodyne and, optionally, 
the defocused homodyne schemes, the light gathered 
from the condensate is assumed to be passed through a 
band-pass filter that suppresses Fourier components out- 
side the interval [ct;iow , ^high] , creating 



#it(r,t) 



1 /*^high r poo 

J i0^c^^T L J —oo 



(C3) 



which is then used to create the interference image (C2). 
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